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INTRODUCTION 


Plasmas  and  fluids  are  rich  in  nonlinear  phenomena,  the  occurrence 
of  instabilities,  the  particle  orbits  in  strong  electromagnetic  field, 
the  weak  and  strong  turbulence,  etc.,  are  just  a  few  familiar  topics  that 
belong  to  this  category.  In  the  past  decade,  many  breakthroughs  occurred 
in  the  understanding  of  nonlinear  phenomena  due  in  large  part  to  the 
advances  in  computing  facilities.  The  discovery  of  solitons  and  the 
subsequent  analytic  solution  method,  the  so-called  inverse  scattering 
method  in  solving  the  nonlinear  wave  equations  completely,  is  one  of  the 
major  events.  Stochastic  dynamics  of  particle  orbits  and  the  onset  of 
turbulence  in  simple  chaotic  dynamical  systems  is  another.  Our  contract 
with  the  Office  of  Naval  Research  carries  out  studies  in  this  field  and 
had  made  many  contributions  to  the  understanding  of  this  exciting  area  of 
research. 


PROJECTS  COMPLETED 


We  summarize  our  accomplishments  under  this  grant  in  this  section. 

It  covers  a  wide  spectrum  of  topics.  However,  we  will  classify  them  into 
main  areas  and  under  each  area  we  will  outline  the  work  done. 


A.  Nonlinear  Internal  Waves  In  Stratified  Fluids 


Internal  waves  that  propagate  in  stratified  fluid  is  a  very 
important  topic  in  physical  oceanography.  It  Is  estimated  that  about 
fifty  percent  of  kinetic  energy  in  the  ocean  existed  in  the  form  of 


internal  waves  and  yet  nonlinear  evolution  of  these  waves  are  not  well 
understood  at  all.  We  have  studied  many  aspects  of  the  nonlinear 
phenomena  associated  with  the  internal  waves.  Topics  Included  soliton 
propagation  in  deep  fluids,  the  Benjamin-Ono  equation,  the  intermediate 
depth  fluid  equation,  the  chaotic  motion  of  solitons,  the  coupling  of 
Internal  waves  with  surface  wa^es  through  a  test  wave  Hamiltonian  model, 
etc. 

1 .  Internal  Wave  Solitons 

It  was  first  derived  by  3enjamin-0no  that  nonlinear  internal  waves 
propagating  in  a  deep  two-fluid  system  is  described  by  a  nonlinear 
integro-dif ferenticl  equation  now  bearing  their  names,  the  Benjamin-Ono 
equation. 

q  +  2qq  +  Hq  =  0  , 
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is  the  Hilbert  transform  of  q.  Benjamin  noticed  the  existence  of  a 
single  solitary  wave  solution 
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and  was  verified  experimentally  by  Davis  and  Acrivos  in  a  water  tank 
experiment.  However,  the  general  solution  of  this  equation  is  not  known 
We  noted  the  numerical  experiment  of  Melss  and  Pereira  that  two  Benjamin 


Ono  solitons  seem  to  maintain  their  identity  after  colliding  with  each 


other,  discovered  a  way  to  construct  general  N-soliton  solutions  of  this 
system.  Later  on,  inverse  scattering  transforms,  infinite  sets  of 
conservation  laws,  Backlund  transformations,  etc.,  were  also  discovered. 
Eventually,  the  complete  solution  is  obtained.  N'ot  only  that,  we  later 
found  that  damping  and  source  can  he  added  into  the  equation  and  yet 
still  permit  N-soliton  solutions.  However,  the  motion  of  the  solitons 
become  chaotic  exhibiting  the  coexistence  of  both  the  coherence  and  chaos 
phenomena. 

When  the  fluid  depth  is  not  infinite,  Ko  et  al.  derived  another 
nonlinear  evolution  equation  to  describe  the  nonlinear  internal  wave. 

qt  +  2qqx+iqxTqxx  =  0  * 


where 


w 

Tq(x)  3  Yg  P  /  coth  — q (y )dy  . 


This  intermediate  wave  equation  is  especially  interesting  because, 
depending  on  the  depth  parameter  6,  it  approaches  to  either  the  KdV  (  g 
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equation  in  the  shallow  water  limit  (5  -*•  0) 
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In  the  infinite  depth  limit  ( 5  +  ®) 


Using  the  Hirota's  method,  we  again  found  the  N-soliton  solution  of 


the  intermediate  equation.  Later  on,  the  inverse  scattering  operator  was 
found  and  the  equation  was  solved  c>mpletely. 

2.  Generation  of  internal  waves  from  its  coupling  to  the  surface  waves 

We  studied  systematically  a  model  system  describing  the  interaction 
of  an  Internal  wave  with  surface  waves  in  the  form  of  multiple-triads. 

In  each  triad,  the  internal  wave  is  a  daughter  wave  and  different  triads 
interact  only  through  the  sheared  internal  wave.  This  test  wave 
Hamiltonian  model  was  first  proposed  and  studied  by  Watson,  Meiss,  and 
Cohen.  They  claim  that  in  the  case  of  two  triads  the  system  is 
integrable,  therefore,  nonstochastic.  One  should  not  expect  the 
therraalization  of  energy  among  the  participating  wave  triads.  However, 
paradoxically,  they  treated  the  system  statistically  by  averaging  over  an 
ensemble  of  initial  data  which  would  then  predict  a  result  of,  for 
example,  rapid  decaying  of  correlation  functions  in  time.  But  in  the 
case  of  a  single  event  we  have  demonstrated  numerically  the  contrary. 

The  correlation  functions,  in  general,  do  not  settle  down  quickly. 

Random  phase  approximation  is  found  to  be  Inadequate.  The  reason  for 
this  is  the  closeness  of  the  system  to  an  integrable  one,  which  retains  a 
long  time  correlation  during  a  large  portion  of  their  trajectory  in  phase 
space. 

B.  Vortices  in  A  Two-Dimensional  Guiding  Center  Plasma 

A  plasma  Immersed  in  a  strong  uniform  magnetic  field  would  be 
constrained  to  move  along  the  magnetic  field  lines  to  form  charged 


rods.  These  charged  rods  tn  turn  would  Interact  with  each  other  and 
result  in  a  drift  motion  which  resembles  the  vortices  motion  in  a  two- 
dimensional  incompressible  fluid.  A  salient  feature  of  these  two- 
dimensional  fluid  systems  is  the  prediction  of  L.  Osanger  of  a  negative 
temperature  equilibrium  state.  Using  a  most  probable  state  approach, 
Joyce  and  Montgomery  derived  an  equation  of  state  of  the  system  which  is 
a  nonlinear  elliptic  partial  different  equation 

2  2 

7  $  +  X  sinh  4=0, 
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where  X  is  positive  for  a  negative  kinetic  temperature  and  <j>  is  the 
electric  potential. 

Due  to  the  nonlinear  nature  of  this  equation,  analytic  solution  is 
hard  to  obtain.  McDonald  et  al.  used  a  generalized  Newton's  method  to 
obtain  several  solutions  of  this  equation  numerically  and  tried  to 
explain  the  formation  of  a  large  vortex  as  the  thermal  equilibrium  state 
We  adopted  instead  the  analytical  method  developed  recently  to  obtain 
periodic  solutions  of  nonlinear  wave  equations  to  this  case  and 
successfully  obtained  the  various  possible  solutions.  Most  importantly, 
the  analytic  methods  allow  us  to  be  sure  of  having  obtained  all  possible 
solutions. 

To  simulate  the  curvature  in  a  real  magnetic  confinement  device,  a 
gravity  is  Introduced  into  the  problem.  The  most  interesting  result 
obtained  is  the  fact  that  in  equilibrium  the  heavier  ion  vortex  would 
float  atop  the  lighter  electron  vortex.  This  exotic  picture  is  due  to 
tiie  negative  temperature  nature  of  the  system.  An  intricate  bifurcation 
diagram  is  also  obtained  which  shows  various  branches  of  solutions 
bifurcate  out  of  and  merge  wttli  each  other. 


.  Parametric  Instabilities  In  a  Plasma 


Parametric  instabilities  were  investigated  in  various  plasma 
conditions. 

1.  Parametric  Instabilities  of  Electron  Cyclotron  Waves  in  Plasmas 


The  recent  developments  of  powerful  microwave  sources  in  the 
millimeter  range  have  caused  substantial  efforts  on  the  electron 
cyclotron  heating  of  large  devices,  e.g.,  Elmo  bumpy  torus,  tokamaks,  and 


mirrors.  We  studied  the  decay  of  an  electron  cyclotron  pump  into  two 
Bernstein  modes  at  the  second  harmonic  cyclotron  layer  in  EBT.  We  find 
it  able  to  account  for  the  heating  of  the  ring  in  the  initial  phase.  The 
coupling  coefficient  for  this  decay  vanishes  for  a  dipole  pump,  whereas 
the  convective  threshold  with  finite  kQ  is  approximately  200  W/cm  .  The 
calculations  could  also  be  useful  in  the  ionospheric  profile  modification 
experiments. 

2.  Cyclotron  Radiation  from  a  Relativistic  Electron  Beam  in  a  Static 


Magnetic  Field 

To  understand  the  nonlinear  stage  of  a  gyrotron,  we  studied  charged 
particle  motions  in  a  self-consistent  electromagnetic  field  and  a 
static  external  magnetic  field.  This  work  is  motivated  by  the  recent 
Interest  on  the  electron  cyclotron  maser  radiation  mechanism  and  the 
advance  in  the  understanding  of  nonlinear  dynamical  systems.  The  full 
nonlinear  dynamics  of  the  particles  are  followed  numerically.  Using  the 
method  of  Poincare's  return  map,  we  found  the  particle's  motion  become 
violent  and  erratic  when  a  critical  energy  is  exceeded,  while  below  this 
threshold  energy,  the  particle  trajectory  is  smooth  and  regular.  The 


transition  from  regular  to  stochastic  behavior  is  very  sharp.  It  may 


have  important  consequences  on  the  limitation  of  masering  devices 


t 


may  also,  on  the  other  hand,  provide  a  way  to  generate  broad  band 
radiations. 

3.  Decay  Instability  of  an  Ion-Cyclotron  Wave  in  a  Plasma 

An  electromagnetic  Ion  cyclotron  wave,  around  the  second  harmonic 
cyclotron  frequency  of  majority  ions  in  the  central  region  of  a  tokamak, 
is  seen  to  decay  efficiently  into  two  Bernstein  waves.  The  coupling  for 
this  decay  vanishes  for  a  dipole  pump;  hence,  a  rigorous  kinetic  theory 
for  decay  instability  of  a  finite  wave  number  pump  is  developed.  For  a 
pump  wave  around  the  cyclotron  frequency  of  majority  ions,  nonlinear 
Landau  damping  and  oscillating  two-stream  instability  turn  out  to  be  very 
important. 

4.  Parametric  Instability  Near  the  lon-Cvc lot  ton  Frequency  in  an 
Inhomogeneous  Plasma 

A  vigorous  kinetic  theory  for  the  parametric  instabilities  of  ion- 
cyclotron,  drift-cyclotron,  and  drift-cyclotron  loss  cone  modes  of  finite 
wave  number  in  an  inhomogeneous  plasma  is  developed.  The  Vlasov  equation 
is  solved  in  the  guiding  center  coordinates  for  two  types  of  equilibria, 
viz.,  Maxwellian  and  loss  cone  velocity  distributions.  For  flute  modes, 
the  coupling  coefficients  are  greatly  simplified.  Ion-cyclotron  waves 
are  seen  to  decay  efficiently  into  low  frequency  drift  waves  in  tokamaks. 
In  a  mirror  machine  a  high  amplitude  drift  cyclotron  loss  cone  mode  of 
negative  energy  decays  into  low  frequency  positive  energy  modes,  giving 


rise  to  explosive  Instability 


In  a  slightly  resistive  magnetized  plasma  Taylor  argued  that 
magnetic  helicity  is  approximately  conserved  while  the  magnetic  energy  of 


the  system  dissipates  away.  The  plasma  would  reach  a  stable  equilibrium 
called  the  minimum  energy  state  which  is  usually  .-re-free.  We  have 
shown  the  existence  of  this  force-free  state  in  lindrical  box.  The 

axial  symmetric  state  which  has  the  minimum  enei  ••’’en  the  aspect  ratio 
L/R  <  1.67  becomes  unstable  to  an  asymmetrical  st.«i.e  as  the  aspect  ratio 
increases  beyond  the  above  critical  value.  This  result  may  have  a 
significant  effect  on  the  construction  of  a  compact  magnetically  confined 
thermonuclear  device  called  the  spheromak. 

2.  Mult i-Sollton  Solutions  of  a  Derivative  Nonlinear  Schrodlnger 
Equation 

In  our  earlier  study  of  lntegrability  of  nonlinear  Hamiltonian 
systems,  a  formalism  was  proposed  to  test  the  lntegrability  of  nonlinear 
wave  equations.  In  particular,  an  example  was  introduced  which  was  shown 
to  be  integrable  by  the  method  of  inverse  scattering.  In  this  paper,  we 
applied  the  Hirota's  method  to  construct  explicitly  the  solutions  of  this 
equation,  a  nonlinear  derivative  Schrodlnger  equation.  Because  of  its 
simple  form,  it  should  find  application  in  many  nonlinear  theories  of 
wave. 

3.  A  Nonlinear  Equation  for  Bidirectional  Ion-Sound  Wave 

Taking  into  account  ion  and  electron  nonlinearity  and  the  effects  of 

finite  ion  temperature,  we  derived  a  nonlinear  equation 

2 

U  -  U  =  U  +  (U  )  +  aU  U  ,  which  can  have  ion-sound  waves 
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propagating  in  both  directions.  Having  a  different  dependence  on  space 
and  time  derivatives,  this  new  equation  describes  more  properly  the 
interaction  of  waves  travelling  in  different  directions.  We  found  two 
soliton  solutions  of  this  equation  from  !ir>ti's  method.  Hut  three  or 
more  than  three  soliton  solutions  are  demonstrated  to  he  nonexistent. 

This  serves  as  an  interesting  physical  example  for  which  only  two  but  no 
more  than  two  solitons  exist. 

4 .  Cylindrical  Ion- Acoustic  Waves 

Nonlinear  ion  acoustic  waves  in  cylindrical  geometry  is  described  by 
a  cylindrical  Korteweg-de  Vries  equation.  Nome r i ca 1 1 y ,  Maxon  and 
Viecelli  showed  that  the  cylindrical  soliton  exists  with  a  till.  These 
solitons  seem  to  be  independent  of  each  other.  They  maintain  their 
identity  after  collisions  with  each  other.  We  showed  analytically  that 
these  solitons  are  indeed  independent  nonlinear  degrees  of  freedom. 
Because  of  the  Airy  functions  that  show  up  in  the  solutions,  tails  always 
exist  which  propagate  in  a  direction  opposite  to  th •  main  part  of  the 
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